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Abstract. Consider a Riemannian metric on two-torus. We prove that 
the question of existence of polynomial first integrals leads naturally to a 
remarkable system of quasi-linear equations which turns out to be a Rich 
system of conservation laws. This reduces the question of integrability 
to the question of existence of smooth (quasi-) periodic solutions for this 
Rich quasi-linear system. 



1. Introduction 

In this paper we study the problem of existence of integrals which are 
homogeneous polynomials with respect to momenta for geodesic flows on 
2-torus. We show that this problem is equivalent to the problem of finding 
periodic solutions for a remarkable quasi-linear system of equations. We 
shall prove that in a domain of hyperbolicity this system can be represented 
in Riemann invariants, and also have a form of conservation laws. Such 
systems are called the rich systems (or semi-hamiltonian) . It is a challenging 
problem to show that a given system of size greater than 2 do not have 
smooth solutions. For Rich systems, the classical analysis by Lax along 
characteristics can be used to establish the shock formation provided the 
so-called genuine non-linearity of eigenvalues is satisfied (see [15] and also [2] 
where this analysis was performed for a problem of integrability for a system 
of 1,5 degrees of freedom). However in our case this genuine nonlinearity 
condition does not necessarily holds. This leaves a hope to find new examples 
of non-trivial integrals for geodesic flows on 2-torus. It would be interesting 
to study if the generalized hodograph method in the region of hyperbolicity 
|16j could give a non-trivial information on smooth solutions for this system. 
We hope to come back to these questions in subsequent paper. 

Consider a geodesic flow on the 2-torus T 2 = M 2 /Z 2 . Let 

2 2 

ds 2 = ^2 gij(q) d q l dq 3 ,H = ^ g lJ (q)PiPj 

i,j=l i,j=l 

be a Riemannian metric and the corresponding Hamiltonian function of the 
geodesic flow. The geodesic flow is called integrable, if the Hamiltonian 
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system 

on T*T 2 ~ T 2 x M 2 , admits a smooth function F(q,p) (called a first integral) 
(functionally independent with H almost everywhere) which has constant 
values along the trajectories of the geodesic flow, i.e. 

We shall denote by r the time variable along the trajectories and -4- will be 
denoted (untraditionally) by dot while the variable t will be used later as a 
coordinate on 2-torus. In this paper we shall deal with the case when F is 
a homogeneous polynomial of certain degree with C 3 coefficients. 

There are two classically known classes of metrics with integrable geodesic 
flows. They are written in a conformal coordinates as follows: 

1. ds 2 = f(x)(dx 2 + dy 2 ), 

2. ds 2 = {f{x)+g{y)){dx 2 + dy 2 ) 

In the first case the Hamiltonian system admits a one parametric group 
of symmetries and has the integral which is a polynomial of the first degree 
with respect to momenta. While in the second case the integral appears to be 
of the second degree and is related to separation of variables in Hamilton- 
Jacobi equation. The existence of the metrics on T 2 with the integrable 
geodesic flows having polynomial integrals of the degree higher than 2 and 
non-reducible to the integrals of the 1 and 2 degree, is not known. Amazingly 
there exist non-trivial examples of geodesic flows on 2-sphere with integrals 
which are homogeneous polynomials of degrees 3 and 4. These examples 
(see [2] and [8]) were inspired by integrable cases of Goryachev-Chaplygin 
and Kovalevskaya in rigid body dynamics (see [6]). 

It is important to mention that in the classical examples 1,2 as well as 
in the examples on the sphere mentioned above the metric is presented in 
conformal coordinates. This approach goes back at least to Darboux, and 
was studied extensively (see for example [11] . [12j . [9]) We the reader refer 
to a nice exposition [15] . 

The main idea of our approach is to work in different coordinates which 
are angle coordinates associated to an invariant torus of the geodesic flow. 
The advantage of these coordinates is a possibility to write the quasi-linear 
equations on the coefficients of the unknown integral and the metric in the 
form of evolution equations as it is described in the following: 

Theorem 1.1. Suppose that the Hamiltonian system (1) has an integral 
F , which is a homogeneous polynomial of degree n. Then on the covering 
plane M? there exist the global coordinates (t,x), where the metric has the 
following form 

ds 2 = g 2 (t,x)dt 2 + dx 2 , 
and the integral F can be written in the form 

k=0 y 
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Where the last two coefficients can be normalized to be a n _i = g and a n = 1. 
Then the commutation relation {F, H} = is equivalent to the system of n 
quasi-linear equations on the unknown U = (ao- 
a n _i = g and a n = I): 

U t + A(U)U x = 0, 
where the matrix A has the form: 



,a n ^2,o,n-i) T (remember 



(2) 



/ o 

O-n-l 





V 











O-n-l 






oi \ 

2a2 — nao 

3d3 — (n — l)a\ 

a-n-i (n - l)a n _i - 3a n _ 3 

a n _i na„ — 2a„_ 2 y 



(3) 



T/te functions a%,g are periodic on the variable x, and quasi-periodic on the 
variable t. 

Let us recall that a hyperbolic diagonal system 

{n) t + A»(ri, . . . ,r n )(n) x = 0,i = l, ..,n 
is called Rich if the eigenvalues satisfy the following conditions 



■a- 



A 7 ; — Ao 



3 r 



g r fc Aj 

A, 



~KJ / V Afe /Vj 

It can be proved that a diagonal system is Rich if and only if it can be 
written in some coordinates as a system of conservation laws (we refer to 
|15j for a nice exposition of rich systems including the blow up analysis and 
to [7].|16j for Hamiltonian formalism for system of Hydrodynamics type). 
For our system the following theorem holds 

Theorem 1.2. The system (2) is the Rich quasi-linear system system. More 
precisely: 

1. In the region of hyperbolicity (all eigenvalues are real and distinct) there 
exists a change of variables (Riemann invariants) (ao, • • • , a n ) — > (r±, . . . , r n ) 
transforming the system to a diagonal form: 

{n)t + Af(ri, . . .,r n )(ri) x = 0,i = 1, ..,n. 

2. There exists a regular change of variables (ao, . . . , a n ) — > (G\, . . . , G n ) 
suc/i iaai Gi, i = 1, 



(Gj(a , 



, n are conservation laws: 
■ ,a n ))t + (Hi{a , . . . ,a n )) a 



0,i = 1, 



, n. 



We shall see that both Riemann invariants and Conservation laws of the 
theorem are ultimately related to the phase portrait of the geodesic flow. 

There are no general methods for proof of existence or non-existence 
of smooth periodic solutions of the system of quasi-linear equations where 
the eigenvalues can collide and become complex. However in [2], [3] non- 
existence of nontrivial periodic solutions for the so-called Benney chain is 
proved. It relies heavily on the property of genuine non-linearity of smallest 
and largest eigenvalues. It seems however that for our system (2) the gen- 
uine non-linearity condition doesn't necessarily holds. This gives a hope to 
find smooth (quasi-) periodic solutions for the system. 
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Let us mention that, starting at t = with a periodic in x initial data 
for the Cauchy problem which is hyperbolic. Then one can "solve" the 
initial value problem and thus to obtain a metric g and the integral F of the 
geodesic flow on a cylinder S 1 ^) x {\t\ < e}. Of course this is local result in 
t and this metric is not necessarily complete. It would be very interesting 
to understand this "catastrophe" in more details. 

In the Section 1 we prove the first theorem. We split the proof of the 
second theorem: in the Sections 2 and 3 we prove respectively, that the 
system (2) can be represented in the Riemann invariants and in the form of 
conservation laws. 
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2. Fermi coordinates. Proof of Theorem 1.1 

The metric of the form ds 2 = g 2 (t,x)dt 2 + dx 2 is of course very well 
known in Riemannian geometry. It can always be introduced locally near 
a given geodesic. In such a case the coordinates (t, x) are called usually 
Fermi coordinates or sometimes equi-distant coordinates. It turns out that 
for integrable geodesic flow these coordinates can be chosen globally. This 
is done with the help of a regular invariant torus of the flow which projects 
diffeomorphically on the base. The existence of such a torus is rather deep 
fact. More precisely, it is proved in [1] (Theorem 1.6) that if the geodesic 
flow on 2-torus has a non-constant polynomial integral, then there exists an 
invariant Lagrangian torus L C T*T 2 , such, that the mapping 

tt\l ■ L - T 2 

is a diffeomorphism, where tt is the canonical projection of the cotangent 
bundle. Moreover this torus lies in a regular level F = f and it consists of 
periodic trajectories. One can assume that / ^ by perturbing a little the 
level if it happened to be zero. Multiplying F by a constant we can achieve 
/ = 1, this we shall assume in the sequel. Let <pi,ip2 denote the "angle" 
part of the action-angle coordinates on L (see []]) . In these coordinates both 
Hamiltonian flows on L, f T of F and h T of H are linearized. Therefore, we 
can assume without loss of generality that the Hamiltonian fields of H, F 
on L equal in these coordinates to: 

sgrad-ff = (0, 1), sgradF = (t^i,^), 

where uj\, UJ2 are some frequencies. Since tt\l '■ L — > T 2 is a diffeomorphism, 
we take <pi,<P2 as coordinates on the base T 2 . We shall keep for them the 
same notations. The following lemma holds 

Lemma 2.1. In the coordinates (pi,y>2 the metric ds 2 has the form 

ds 2 = a((f\, ip2)d(p 2 + 2bd(fidf2 + d<P2i 
where b € M is a constant and a((pi, 992) is a positive periodic function. 
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Proof. Let us assume that in the coordinates ip\ , ip2 the metric ds 2 has the 
form 

ds 2 = a(ipi,(f 2 )d(pf + 2b((pi,(p 2 )dipid(p2 + c(<pi,<p2)d<p%. 
According to our choice of coordinates ((pi,ip 2 ) the curves, defined in the 
parametric form (tpi = (fio, ip 2 = 9220 + r )> are geodesies. Since the length of 
a velocity vector of these geodesies is one, we get 0(^1,^2) = 1- Moreover, 
we have 

< 7r*(sgradif), 7T* (sgradF) >= bu)\ +u>2, 

on the other hand it can be easily checked using the homogenuity of F, that 
on L the following holds: 

< 7T* (sgrad-ff ) , 7r* (sgradF) >= pdq(sgr&dF) = 

dF dF 

= Pit; \~ P2^~ = nF = nf = n. 

dpi OP2 

Here pdq stands for the canonical one form in the cotangent bundle. Thus 
bu>i + UJ2 = nf = n, and consequently, b = (n — a^V^i is a constant. 

□ 

Let's introduce new coordinates (t,x): 

(pi = t, ip 2 = x — bt. 
Therefore, in the new coordinates the metric takes the form 

ds 2 = g 2 dt 2 + dx 2 , 

where g 2 = a—b 2 . In the coordinates t, x the vector field 7r* (sgradi 7 ) becomes 
equal to (u>i,uj2 + bu>i) = 

Write the integral F with respect to these coordinates: 

z? — a k(t,x) k k 

t ~ Z> Q n-k Pi ^2, 
k=0 y 

where pi,P2 are the conjugate momenta to t, x respectively. Then obviously 
for the torus L we have L = {p\ = 0,p2 = 1} and so F\l = a n = 1. 
Moreover, partial derivatives of F are easily computed on L to be: 



dF 

dpi 



n-l dF_ 
9 ' dp 2 



na„ = n. 



Compare these values with ^(sgradF) = (a>i,n), we conclude 

= UJi. 

9 

In addition, we shall assume that u± = 1. Indeed, this can be achieved easily 
by a rescaling the t variable, t — > rt for some constant r. In order to complete 
the proof of Theorem 1 1 . 1 1 one has to write the condition {F, H} = explicitly 
and in this way to get system (2). We omit this computation. Notice that by 
the very construction of the coordinates t, x it follows that for any function 
/ on the torus T 2 it can be written as a periodic function f{<pi,<p2), which 
equals f(t, x — bt). Therefore, / becomes periodic in x and (-quasi) periodic 
in t. Theorem 1.1 is proved. 
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3. RlEMANN INVARIANTS 

In this Section we show that the system (2) can be represented in Riemann 
invariants in the hyperbolic domain by change of field variables (ao, • • • , a n ) — > 
(n, ■ ■ ■ ,r n ). 

2 

Let's fix the energy level of the Hamiltonian H = = \- Assume 

that p\ = gcosip, p2 = simp, where <p is an angular coordinate in the fibre. 
Then F = F(t, x, ip) becomes a trigonometric polynomial. We have 

n 

F = F(t,x,<p) = y^ j a k cos n ~ k ip sin k ip; a n -i = g,a n = l. (4) 

fc=0 

dF . cos ifi , , 

— = F t t + F x x + F^ip = F t + F x simp + F^ip = 0. (5) 

dr 9 

Lemma 3.1. The following equality holds 

cos™ tp 

where xa is the characteristic polynomial of the matrix A with the relation 
A = g tan ip. 

The proof can be easily obtained by a direct verification. Remarkably 
this lemma gives a very clear geometric meaning of the eigenvalues of the 
matrix A(U) of our system. They correspond precisely to those ip where the 
derivative i 7 ^, i.e where the invariant tori of the geodesic flow are tangent 
to the fibres. 

With the help of the lemma and the equations (4), (5) one can easily write 
the system (2) in the form of Riemann invariants. Let Aj, i = 1, . . . , n be n 
distinct real eigenvalues of A(U). Define 

7T 7T 

ri = F((pi), <pi = arctan(Xi/g), <£»€[--;-] for i = l,...,n. 

It follows from the lemma that F' (ipi) vanishes and so by the equation (5) 
it follows that 

(n)t + \{ri) x = 0, % = 1, ..,n. 

Moreover formula (4) implies that (n, . . . , r n ) are regular coordinates, since 
the Jacobian matrix for this change of variables equals 

^ Tj-^J = (cos n ~ k ipi sin k ip^j , i = 1, . . . , n; k = 0, . . . , n — 1 

which is non-degenerate matrix if ipi £ [ — f ; § ] are distinct. 



4. Conservation laws 

In this Section we show that the system (2) can be represented in the form 
of the conservation laws. The method of proof of this statement is based on 
the work [5] . Let us first rewrite Hamiltonian system of the geodesic flow as 
a system with 1,5 degrees of freedom away from the invariant torus L. This 
can be done as follows. Any geodesic 7(r) different from the "vertical" ones 
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{t = const} has to have ^ 7^ and so can be written as a graph {x = x(t)}. 
The length functional can be rewritten in the following way: 



sd) = J\i\dr=J ^ 2 (f) 2 +(^)V = J v^rr? = J lm, 



1 1 i ill Hi u ji li 1 1 L \ ; .1 

4^. Then the momenta variable for x' is p = 44 



where £(x',x,t) = \J (x') 2 + g 2 where we write (untraditionally) x' = ^jf , 7 



Then Legendre transform for £ is x, t) = px' — C = —g\/l — p 2 . So 
Hamiltonian equations for geodesic flow get the form: 



, dU gp 

x = 



dp y/T^ 



&H 



P = —q^ = 9xV 1 -P ■ 

Notice that in coordinates (p, x, t) formula (4) of the previous section for F 
gets the form 

n 

F(p,x,t) = Y,*k(l-p 2 )-*p n - k , (4') 

k=0 

where the dependence on (t, x) enters through the coefficients only. 

Assume L\ is any torus which is invariant under the flow given as a graph 
of a function p = f(t,x). Then the invariance condition means that the 
1-form f t dt + f x dx — dp vanishes on the Hamiltonian vector field sgradW. 
Therefore we have: 

f + f 8H dH_ 
Jt + Jx-z. — I- ^— — u, 
op ox 

which is equivalent to 

% + -^n(f(t,x),x,t)=0. (6) 

Equation (6) implies that every invariant torus which projects diffeomor- 
phically to the base produces in fact a conservation low for our quasi-linear 
system (2). More precisely, let us choose n disjoint tori L±,... ,L n in the 
neighborhood of the torus L, lying in regular levels {F = q} such that their 
projections to T 2 are diffeomorphisms. 

Each of these torii Lj is a graph of a function fi, p = fi(ao, ■ ■ ■ ,a„_i). 
They satisfy the equation 

F(ao,...,a n -i,p) = Ci. (7) 

Notice that fi depend on (t, x) implicitly through the coefficients. Using 
(6) for fi we get the conservation laws: 

f + !;(«(/.,*,*)) = o. 

Let's show that the transition from the variables ao, . . . , a n _i to the variables 
/i, ...,/„ is a diffeomorphism. Introduce for convenience 



ui = a , • • • , u n -i = a n -2,u n = a n _i = g. 



MISHA BIALY AND ANDREY E. MIRONOV 



We have 
Therefore 



F(m, . . .,u n ,fi) = c 



OF dF 

duj dp 



P =fi duj 

or equivalently in the matrix form we have: 



= 



where 



dF it 
dm Ui 



A + B 



9h 

dfn 



dfi \ 

dfn 

du„ ) 



dF 

du„ 



(fn) 



/ 9F 



\ 



P=fi 





dF 
dp 



P=fn 



) 



Note that the matrix A = (^§£:(fi)j 1S essentially Vandermonde matrix, 
in particular it is non-degenerate since the torii Lj are disjoint. The diagonal 
matrix B in the formula is non-degenerate either, because the torii Lj were 
chosen to be regular graphs. Therefore, 

d(fi, ■■■ ,f n ) _^ 
d{ui,...,u n ) 

This completes the proof of Theorem 2. 

In practice it's difficult to reverse the equation (7) explicitly. Therefore, 
instead of solving of the equation (7) we act in the following way. Torus L 
has the the equation p = 1 in the energy level. The value of the integral F 
equals 1 on L. For a small parameter e solve the equation 



P )2p 



n—k 



l + e 



k=0 



by substituting the power series 

p = 1 - G 2 e 2 - G 3 e 3 .... 

Then it follows from (6) that all the coefficients Gj(ao, • • • , a n -i) are con- 
servation laws. 

Example 1. Let n = 3. In this case U = (ao, ai, a2) T , with a 2 = g,a 3 = 1. 
The matrix A{U) has the form: 

ai 
A(U) = [ a 2 2a 2 - 3a 
a2 3a3 — 2ai 



The integral has the form: 

F = oo(l -P 2 )^ +oi(l -p 2 )p + a 2 (l -p 2 )^ +a 3 p 3 . 

Write 

p = 1 - G 2 e 2 - G 3 e 3 - G 4 e 4 . . . 



(8) 
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and substitute into the expression for F in (8), and equate to 1 + e. Com- 
puting coefficients for e, e 2 ,e 3 and equating them to f,0,0 respectively one 
gets the following conservation laws: 



One can prove in general for any n, that all the conservation laws obtained 
in this way are rational functions of the field variables. It is not an easy 
exercise to verify by hands that these are in fact conservation laws. 



1. An intersting direction for further study would be to understand alge- 
braic properties of the system. Notice that the matrix A(U) depends linearly 
on the components of U, this case was studied first by Gelfand, Dorfman 
see HQ]. 

2. It is important to understand what the analysis along characteristics 
can give for smooth solutions of the system. 

3. It is not clear to us if the generalized hodograph method can give a 
nontrivial information for system (2). 

4. Let us mention an interesting connection between real eigenvalues of 
the matrix A(U) and the so called separatrix chains introduced in [3]. The 
geometric language seems to be of tight relation with analytic properties of 
our system. 
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